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We theoretically study measurement induced-dephasing of a superconducting qubit in the circuit 
QED architecture and compare the results to those obtained experimentally by Schuster et al, 
[Phys. Rev. Lett. 94, 123602 (2005)]. Strong coupling of the qubit to the resonator leads to a 
significant ac-Stark shift of the qubit transition frequency. As a result, quantum fiuctuations in the 
photon number populating the resonator cause dephasing of the qubit. We find good agreement 
between the predicted line shape of the qubit spectrum and the experimental results. Furthermore, 
in the strong dispersive limit, where the Stark shift per photon is large compared to the cavity decay 
rate and the qubit linewidth, we predict that the qubit spectrum will be split into multiple peaks, 
with each peak corresponding to a different number of photons in the cavity. 
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I. INTRODUCTION 



Superconducting qubits are promising building blocks 
for the realization of a quantum computer Sev- 
eral experiments have shown coherent control of a sin- 
gle qubit 0, 1MB; 1MB and two-qubit experiments have 
been realized fn. IMIqIIioI ITU . Recently, it was suggested 
that superconducting qubits can be strongly coupled to 
distributed or discrete LC circuits in a way that opens 
the possibility to study quantum optics related phenom- 
ena in solid-state devices [H El 111 El Eg. This 
concept has been successfully demonstrated experimen- 
tally 0, El El ES is and effects associated with the 
quantum nature of the microwave electromagnetic field 
have now been seen in the form of vacuum Rabi split- 
ting El and measurement induced dephasing via pho- 
ton shot noise E3- In this paper we present a detailed 
analysis of the quantum fluctuations of the photon num- 
ber in the cavity and its effect on the qubit spectrum. 
We also show that access to the extreme limit of strong 
dispersive coupling should allow direct observation of the 
photon number distribution in the cavity. 

An advantage of some of these circuit QED analogs 
of cavity QED is that the cavity presents a well defined 
electromagnetic environment to the qubit which can lead 
to enhanced coherence times of the qubit 0| . This well 
defined environment makes quantitative predictions for 
superconducting qubits more tractable. This was shown 
in Ref. 01 where we have studied Rabi oscillations in a su- 
perconducting qubit strongly coupled to a superconduct- 
ing transmission line resonator. Due to the detailed un- 
derstanding of the measurement process, we were able to 
make quantitative predictions about the measured pop- 
ulations in the Rabi oscillations and observe high visibil- 
ity fringes. Moreover, as we showed in Ref. |l9f. pop- 



ulating the strongly coupled resonator with a coherent 
microwave field can lead to a significant ac-Stark shift of 
the qubit, even in the situation where detuning between 
the cavity and qubit frequencies is large. Due to the 
shot noise in the number of photons populating the res- 
onator, this ac-Stark shift leads to measurement-induced 
dephasing of the qubit. This is similar in spirit to the 
experiment on Rydberg atoms in a 3D cavity reported in 
Refs. |2l|, 123, m |2J| . In those time-domain experiments, 
the visibility of Ramsey fringes was shown to decay with 
an increase of the strength of dispersive coupling to the 
cavity. 

In this paper, we will expand on the theoretical model 
presented in our experimental paper |l9l| (hereafter re- 
ferred to as the Letter) where we observed the AC stark 
shift and measurement induced dephasing in a circuit 
QED device. We will start in Sec.mwith a brief review of 
the important features of circuit QED. In Scc. lIIII thc ex- 
perimental results reported in the Letter will be reviewed. 
We then present in Sec.lVltwo theoretical models describ- 
ing measurement-induced dephasing. We first start with 
a simple model which assumes Gaussian fiuctuations of 
the qubit's phase. This is the model that was briefly pre- 
sented in the Letter to explain the experimental results. 
We then present a more gen eral approach based on the 
positive P-representation |25( which does not require the 
Gaussian assumption. For the experimental parameters 
reported in the Letter, these two approaches give identi- 
cal results. However, in the limit of strong coupling and 
very high Q resonators, the second approach shows that 
qubit spectrum will exhibit structure at several distinct 
frequencies due to the underlying discrete energy levels 
of the cavity. That is, we predict that the qubit spectrum 
will split into multiple peaks, with each peak correspond- 
ing to a different number of photons in the cavity. We will 
refer to this as number splitting of the qubit spectrum. 
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Experimental observation of this effect would be a direct 
demonstration of number quantization in the dispersive 
regime. We also discuss how, by using irradiation which 
is off-resonant from both the cavity and the qubit, one 
can obtain substantial ac-Stark shifts without significant 
dephasing and how this could be used as the basis of a 
phase gate for quantum computation. 



II. CAVITY QED WITH SUPERCONDUCTING 
CIRCUITS 

A. Jaynes-Cummings interaction 

In this section, we briefly review the circuit QED ar- 
chitecture first introduced in Ref. |l3l | and experimentally 
studied in Refs. shown in Fig. Hthe sys- 

tem consists of a superconducting charge qubit Ill [l^ l2q 
strongly coupled to a transmission line resonator |27l |. 
Near its resonance frequency w^, the transmission line 
resonator can be modeled as a simple harmonic oscilla- 
tor composed of the parallel combination of an inductor 
L and a capacitor C . Introducing the annihilation (cre- 
ation) operator d*-^', the resonator can be described by 
the Hamiltonian 



(2.1) 



with u!,. = \j\fLC. Using this simple model, one 
finds that the voltage across the LC circuit (or, equiv- 
alently, on the center conductor of the resonator) is 

Vlc = Kms(^^ + a)i where V^^g = y/hu\j2C is the rms 
value of the voltage in the ground state. An important 
advantage of this architecture is the extremely small sep- 
aration 6^5 /im between the center conductor of the 
resonator and its ground planes. This leads to a large 
rms value of the electric field E°^^ = V^^^Jb - 0.2 V/m 
for typical realizations 01^3 ^3- As illustrated in Fig.^ 
by placing the qubit at an antinode of the voltage, it will 
strongly interact with the resonator through the large 
electric field E^^^. 

In the two-state approximation, the Hamiltonian of the 
qubit takes the form 



E, 



(2.2) 



where E^i = 4£'c(l — 2ng) is the electrostatic energy and 
iJj = iJj^niax cos(7r$/$o) the Josephson energy. Here, 
Ec — e^/2Cs is the charging energy with Cs the total 
box capacitance and Ug = CgVg/2e the dimensionless 
gate charge with Cg the gate capacitance and Vg the gate 
voltage. -Ej,max is the maximum Josephson energy and 
$ the externally applied flux, with $0 the flux quantum. 

Due to capacitive coupling with the center conductor, 
the gate voltage Vg = Vg'^ + Vlc has a dc contribution 
V^'^'^ (coming from a dc bias applied to the input port of 
the resonator) and a quantum part Vlc • When working 




FIG. 1: (Color online) Schematic layout and lumped element 
version of the circuit QED implementation. A superconduct- 
ing charge qubit (green) is fabricated inside a superconducting 
ID transmission line resonator (blue). 



at the charge degeneracy point rig =1/2 where dephas- 
ing is minimized [j| and neglecting fast oscillating terms, 
the resonator pl us q ubit Hamiltonian takes the Jaynes- 
Cummings form [iSj 



fig (d^CT- + (7 -I 



(2.3) 



where Wa = E^/Ti is the qubit transition frequency and 
g = e{Cg/Cs)V^^^/h is the coupling strength. 

As shown in Ref. the qubit can be measured 

and coherently controlled by applying microwaves, of 
frequency uj,-i and Us respectively, to the input port of 
the resonator. This can be described by the additional 
Hamiltonian 



Hn = 



E 

j=s,rf 



TlEj {t) 



■ de+*"^*) 



(2.4) 



where Sj{t) is the amplitude of the external drives at rf 
and spectroscopy frequencies. 



B. Dispersive regime 

In the situation where the qubit is strongly detuned 
from the cavity, |A| = [ajr — LOal 2> g, the total Hamil- 
tonian Hjc + ^^D can be approximately diagonalized to 
second order in g/A to yield the following quantized ver- 
sion of the dynamical Stark shift Hamiltonian [Tsl 



i/eft = TiUrO^a + - [ua + 2xd^d) dz 



j=s,rf 



(2.5) 



j=s,rf 



Here Coa = Wa + X is the Lamb shifted qubit frequency 
and we have defined x = g^ j A. The term proportional 
to a^aa^ can be interpreted as a shift of the qubit tran- 
sition frequency depending on the photon number in the 
resonator (ac-Stark shift) or as a pull on the resonator 
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frequency by the qubit. As will be shown later, quantum 
noise in the photon number a leads to dephasing of the 
qubit. 

Dephasing due to coupling to the cavity field was also 
studied using Rydberg atoms coupled to a 3D microwave 
cavity |2lL l22l |23|. In this experiment, atoms were sent 
one at a time through the cavity and interacted for a 
finite time with the field. The state of the atoms was 
finally read out by ionization . Visibility of the Ram- 
sey fringes was measured as a function of the detuning 
from the cavity, hence as a function of x- Dephasing was 
shown to increase with the strength of dispersive cou- 
pling X to the cavity In this paper, we will instead 
consider dephasing of the qubit due to the resonator field 
by looking at the qubit spectrum as measured by trans- 
mission of the cavity field. 

We note that, in practice, lOs is chosen to be close to 
Cja and the last term of Eq. 1)2. 5|l with ujs causes Rabi 
flopping of the qubit. Moreover, as further discussed be- 
low, we choose Wj-f ~ ujj. — to measure the state of 
the qubit, where is the detuning of the measurement 
probe from the bare cavity frequency. In this situation, 
the last term of Eq. I|2.5|l with oj^f is largely detuned from 
the qubit and does not lead to qubit transitions. As first 
noted in the original proposal by Brune et al. [2^l30| this 
measurement Hamiltonian is therefore highly quantum 
non-demolition |3ll | with respect to measurement of the 
qubit state. Conversely, if the dispersive coupling term 
is dominant in the Hamiltonian, then QND measurement 
of photon number is also possible. Physical implemen- 
tation of QND readout for superconductin g qu bits has 
been achieved in the microwave regime |3. llTl IT^ . 
but values of the coupling % large enough to allow QND 
readout in the dispersive regime for the photon number 
have not yet been achieved in any system. However re- 
markable experiments on Rydberg atoms have achieved 
photon number readout in a non-dispersive (i.e. degen- 
erate) regime |2lll3^ . 



C. Damping 

Coupling to additional uncontrolled bath degrees of 
freedom leads to energy relaxation and dephasing in the 
system. Integrating out these degrees of freedom leaves 
the qubit plus cavity system in a mixed state p{t) whose 
evolution can be described by the master equation |3 ij 

^-^[H,p] + KV[a]p + jMa^]p + ^Via,]p, ^'^'^^ 

where T>[L]p = (^LpD - DLp - pDl^ /2 describes 

the effect of the baths on the system in the Markov ap- 
proximation. The last three terms of Eq. (|2.6|) correspond 
to loss of photons at rate k, energy relaxation in the qubit 
at rate 71 and pure dephasing of the qubit at rate 7^. 



In the dispersive regime, the operators describing en- 
ergy relaxation and dephasing should be transformed in 
the same way as was done in Eq. (|2.5|l . This leads to 
small corrections, of order (g/A)'^, to the master equa- 
tion that are omitted here. 



III. EXPERIMENTAL RESULTS 
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FIG. 2: (Color online) Measured spectroscopic lines (blue 
lines) at (a) intra resonator photon number n ~ 1 with fit 
to lorentzian line shape (solid line) and at (b) n ~ 40 with 
fit to gaussian line shape (solid line). Dashed lines are best 
fits to (a) gaussian or (b) lorentzian line shapes, respectively. 
The qubit transition frequency Qa at low measurement power, 
the half width half max Suhwhm and the ac-Stark shift ojac 
of the lines are indicated. 



In this section, we briefly review some of the exper- 
imental results already presented in the Letter. Only 
the results that are directly discussed in the present pa- 
per will be presented and the details of the experiment 
can be found in Refs. H IH El H^. In the Letter, 
we reported spectroscopic measurements of the qubit as 
a function of measurement power. The qubit spectro- 
scopic line is shown in Fig. |2] for two average photon 
numbers n in the resonator, corresponding to two in- 
put measurement powers. The relevant experimental 
parameters are A/27r = 105 MHz, g/27r ^ 5.8 MHz, 
Ar/27r = 0, and K/27r = 0.57 MHz and a dephasing time 
(in the absence of power broadening) of T2 > 200 ns. 
These values correspond to a relatively small cavity pull 
of x/27r w 0.32 MHz, or x/k ~ 0.56 in units of the cavity 
line width. It is important to note that at these relatively 
low detunings A, there can be a qubit contribution to the 
cavity line width. 

As discussed in the Letter, at low measurement power, 
the line shape of the qubit spectrum is Lorentzian but 
as the measurement power increases, the line shape ap- 
proaches a Gaussian. As shown in Fig.|21 this is also seen 
in the dependence of the half-width at half-maximum 
'^I'HWHM of the qubit line shape on n which goes from 
cx n to cx \ffi as measurement power increases. This fig- 
ure also shows theoretical results that will be discussed 
below. 
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FIG. 3: (Color online) Full, red curve: Measurement broad- 
ened qubit line width (5j/hwhm as a function of the input mea- 
surement power or average photon number as predicted by the 
lowest order dispersive approximation. Green, dotted curve: 
Same as red but taking into account the non-linear reduction 
in the cavity pull and plotted as a function of input power. 
The symbols are the experimental results. Symbols and color 
scheme are described in the text. The vertical line indicates 
the critical photon number ricrit- The parameters are those 
given in section 11111 The blue dashed line is the calculated 
HWHM for A,./27r = 32 MHz. It clearly shows that measure- 
ment induced dephasing is small at large Ar where informa- 
tion about the state of the qubit in the transmitted signal is 
also small. This can thus be used as the basis of a phase gate. 



In the Letter, we have already provided a theoretical 
explanation for this behavior. In section we review 
and expand on this model. We then explain how in the 
limit of very large cavity pull, the results can be signifi- 
cantly different. 



IV. AC-STARK SHIFT 

In the lowest order dispersive approximation [Eq. H2.5|l ] 
the predicted ac-Stark shift o^ac = 2xn will be a linear 
function of the mean photon number, n. However, this 
approximation only holds at low photon numbers and 
breaks down on a sc ale g iven by the critical photon num- 
ber ncrit = This is illustrated in Fig. Hi) 
where the ac-Stark shift, calculated from the lowest order 
dispersive approximation (red solid line) and the exact 
eigenvalues (blue dashed line) of the Jaynes-Cummings 
model ^3 are plotted for the experimental parameters. 
Also shown in this figure is ncrit (vertical blue line) which 
for the experimental parameters is about 82 photons. 
Here we see that as n increases, the exact Stark shift 
begins to fall below the lowest order dispersive approxi- 
mation even before n reaches ricrit- 

In Fig. Qd) the experimental results (soHd blue points) 
are plotted as a function of probe power, Prf (extend- 
ing up to powers larger than those presented in the Let- 
ter). To convert between n and Prf, we assume that 
-Prf — ATiwrfP where p is the photon flux at the res- 



onator and A is a scaling factor that takes into account 
the large attenuation that is placed between the probe 
generator and the resonator (to eliminate black body ra- 
diation). From the lowest order dispersive approxima- 
tion, the average photon number when driving the cavity 
at = is 



(k/2)2+x2' 



(4.1) 



By using the lowest order dispersive approximation for 
the ac-Stark shift and the line of best fit to the experi- 
mental points (in the linear regime at low power) A can 
be determined. Doing this gives the red solid line in 
Figs.0lD) and d). The calibration shows that ricrit occurs 
at w 110/iW. The experimental results clearly show the 
breakdown of the lowest order dispersive approximation. 
The data points fall below the linear prediction but not 
nearly as much as the blue dashed curve in (a) predicts. 
In fact, the data points follow fairly closely the linear in n 
dependence of the lowest order dispersive approximation 
in Eq. (|2.5(l for larger powers than expected (up to and 
well above ricrit)- 

It is possible to understand why the experiment agrees 
with the simple dispersive approximation for larger probe 
powers than expected by considering the following sim- 
ple model. We assume that, at these large powers, the 
ac-Stark shift is still given by the dispersive approxima- 
tion Ixn^ but wc now take into account the non-linear 
cavity pull (which is x a-t low h). From the eigenval- 
ues of the Jaynes-Cummings Hamiltonian the cavity pull 
can be calculated as a function of n. This is shown in 
Fig. H). From this figure, we see that the cavity pull 
reduces as the number of photons in the resonator is in- 
creased. We thus replace x by x(n). The second aspect 
of our simple model is the non-linear dependence of the 
average photon number n with input power P due to the 
power-dependence of the cavity frequency. To account 
for this we simply replace x Eq. (|4.1|) with x(^) and n 
becomes a non-linear function of input power. This non- 
linear dependence of the photon number on the input 
power is illustrated in Fig. 0Jl) as the green dotted line. 
This is a precursor to bistability in this system • Us- 
ing these two expressions, we have for our simple model 
of the non-linear ac-Stark shift 2xiP^n[P^. This expres- 
sion is plotted in Fig. 2b) (green dotted line) with a new 
scaling factor A' « 0.905A calculated by the best fit for 
the experimental data (here we use the complete data 
set). This simple model produces a result that is linear 
for a larger range of powers and is closely consistent with 
the experimental results. It happens that for the partic- 
ular experimental parameters, the two non-linear effects 
almost cancel each other out and result in the green dot- 
ted line being more linear than expected. 

We emphasize that in the Letter, only the low power 
(below ricrit) part of the ac-Stark shift was studied and 
was fit only with the linear dispersive model. Compari- 
son with the results of the non-linear model shown here in 
Fig. ^) shows that the calibration of the cavity photon 
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FIG. 4: (Color online) (a) Full, red curve: ac-Stark shift as a function of the average intra-cavity photon number using the 
lowest order dispersive approximation. Dashed, blue curve: ac-Stark shift as a function of photon number calculated from the 
exact eigenvalues of the Jaynes-Cummings Hamiltonian. (b) Blue dots: experimentally measured ac-Stark shift as a function 
of external microwave input power. Full, red curve: Predicted ac-Stark shift within the lowest order dispersive approximation. 
The conversion factor from photon number to external microwave drive power was determined by fitting the red curve to the 
linear portion of the data at low power. Green, dotted curve: Same as red but taking into account the non-linear reduction in 
the cavity pull (see part c) and the non-linear increase in the average photon number (see part d) with microwave drive power. 
For the particular experimental parameters these two effects almost cancel each other out and result in the green dotted line 
being nearly linear out to much greater input powers than expected, (c) Cavity pull as a function of average photon number 
n. The red solid line is the result of the dispersive approximation (±x) while the dashed blue curve is obtained from the exact 
eigenvalues of the Jaynes-Cummings model, (d) Average photon number as a function of input power. The full red line is the 
result of the lowest order dispersive dispersive approximation Eq. (14.111 fit to the data in (b) at low power. The dotted green 
line is the non-linear model with x replaced by x(^) in El- J4.1II The vertical line in all plots indicates the critical photon 
number ricrit = A^/4(j^ which indicates the scale at which the lowest order dispersive approximation breaks down. For the 
experimental parameters (given in section ITTTl ricrit ^ 82 which corresponds (within the lowest order dispersive approximation) 
to Prf ~ llOAtW. 



number in terms of the drive power is low by approxi- 
mately 50% at the highest power shown in Fig. 5 of the 
Letter. We also emphasize that our treatment here of the 
non-linearities is only approximate. 

V. MEASUREMENT-INDUCED DEPHASING 

By monitoring the transmission of the cavity using 
heterodyne detection, one has access to the average of 
the cavity field (a(t)). As shown in Ref. the phase 
(j){t) = arg{(a(t))} is directly related to the population 



of the qubit ((T^(t)). As a result, by recording the phase 
(^(t) as a function of the excitation frequency tOsi one has 
access to the absorption spectrum of the qubit [23 

1 f°° 

S{oj)^—j dte^-'{a.it)a+ms, (5.1) 

where the subscript s implies that the expectation value 
is taken in the steady state. The dephasing rate can be 
determined through the half width at half maximum of 

Using the quantum regression formula |25l |. the corre- 
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lation function ((T_(t)(T_|_(0))s can be evaluated as 



(a_(t)<T+(0)), 



t > 
t < 0, 



(5.2) 



where is the steady-state density matrix. This allows 
us to rewrite the spectrum as 



Situ) 



1_ 



-Re 



dt e"^'{a^{t)&+{0))s 



(5.3) 



We start by calculating this spectrum by assuming 
gaussian statistics for the qubit's phase noise. This ap- 
proximation is sufficient for the experimental parameters 
considered above but breaks down in the situation where 
the cavity pull is large x ^ Moreover, this simple 
approach has the advantage of presenting the essential 
physics in a transparent way. We then show how to go 
beyond the gaussian approximation by using the positive- 
P hmction approach |25j. 



A. Gaussian approximation for the phase 

As mentioned above, quantum fluctuations Sn in the 
photon number around its average value n will lead to 
accumulation of a random relative phase between the 
amplitudes of the two basis states of the qubit and 
hence to dephasing. We first consider the situation 
where the qubit is prepared in a superposition \ip{0)) = 
{\g) + |e)) /^/2 of its basis states and that the cavity is 
populated by a coherent state with average photon num- 
ber n. As it evolves under the dispersive Hamiltonian 
Eq. 12.5|l . the qubit superposition picks up a relative 
phase factor 



ip{t) = + 2x / dt'n{t'), 
Jo 



(5.4) 



where Wa = + X is the Lamb shifted qubit transition 
frequency. It is convenient to express the second term as 
its mean value plus fluctuations about the mean 

(p{t) = (p + 6ip{t) = Uat + 2xnt + 2x / dt'5n{t'), (5.5) 

Jo 

with n the average photon number in the cavity leading 
to the ac-Stark shift (see Fig. ^ and 5n{t) the random 
excursions about this mean. 

In a frame rotating at the Lamb and ac-Stark shifted 
qubit transition frequency, we obtain for the correlation 
function (t > 0) 



(a_(t)rT+(0))3 = Tr[a_e^*<7+p,] 

= Tr[a_e^*(|e)(g|)] 



(5.6) 



where 72 = This is the off-diagonal component 

of the reduced qubit density matrix. Assuming gaussian 



statistics for the phase Sip(t), the cumulant expansion is 
exact and we obtain 133 



(a_(0^+(0)). 



e '^^^ exp 



-2x' 1 1 dtidt2{5n{ti)Sn{t2)) 



(5.7) 



This expression involves the photon-photon time corre- 
lator which for a two-sided symmetrically damped driven 
cavity takes the form [T3l | 



{Sn{ti)Sn{t2)) ~ ne 



-t2\ 



(5.8) 



leading to 



(a_(t)a+(0)) 
= e"'''^* exp 



2 



1 + exp 



K\t\ 

2 



(5.9) 



where Oq = tan"-"^ 2x/ k « 2x/ k is the magnitude of the 
accumulated phase shift for the transmitted photons due 
to the coupling with the qubit in the small pull approxi- 
mation (x ^ ^) at Wrf = i^r- 

We now consider two simple limits of the above re- 
sult. First, in the situation where the mean cavity pho- 
ton number n is small, fluctuations of the photon number 
will only weakly contribute to dephasing. In this situa- 
tion phase decay occurs on a long time scale with respect 
to 1/k. In this limit, exp(-K|t|/2) w and Eq. H5.9|) 
reduces to 

(<7_(t)CT+(0)),, « exp [- {72 + 2nfiel} \t\] . (5.10) 

In this situation, the measurement induced-dephasing 
only adds to the intrinsic dephasing rate 72. This is 
because, in this long time limit, the phase undergoes a 
random walk process leading to an exponential decay of 
the coherence. The Fourier transform of this expression 
leads to a Lorentzian spectrum with half-width at half 
maximum 72 -l- Fm, where Fm = 2Kn0g is the measure- 
ment induced dephasing rate in the small pull limit (see 
next section). 

On the other hand, in the large n limit, phase decay 
can occur on a time scale much shorter than the cav- 
ity lifetime, t ^ k~^. Indeed, since n is large, a small 
fraction of n leaking out of the cavity in a time t <C 
conveys enough information to infer the state of the qubit 
and hence to dephase it completely [ssj . In this situation, 
expanding exp(— «;|t|/2) in Eq. H5.9|l . we obtain 



(<7_(i)a+(0)), « exp [-72|t| - 2nxh^] 



(5.11) 



The large n limit does not lead to an exponential decay 
and the spectrum will be a convolution of a Lorentzian 
and a Gaussian. This corresponds to inhomogcneous 
broadening of the qubit due to the Poisson statistics of 
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the coherent state populating the cavity. In this situa- 
tion, the half width at half maximum therefore scales as 

Using the full expression Eq. (|5.9|l and moving back to 
the lab frame, we obtain for the spectrum of the qubit in 
the Gaussian approximation for the phase 



^ K 



2^ 3 



"27r ^ 
j 



2^ J 



(5.12) 



where f ^ = 2{'-^2+^m.)+j>i- The spectroscopic line shape 
is given by a sum of Lorentzians, all centered on the ac- 
Stark shifted qubit transition but of different widths and 
weights. 

As expected from the above discussion, we see from 
Eq. H5.12|) that if the measurement rate Fm is much 
smaller than the cavity decay rate k/2, then only a 
few terms in the sum contribute and the spectrum is 
Lorentzian. On the other hand, when the measurement 
rate is fast compared to the cavity damping, the spec- 
trum will be a sum of many Lorentzians, resulting in a 
gaussian profile. In this situation, dephasing occurs be- 
fore the cavity has had time to significantly change its 
state, leading to inhomogencous broadening as discussed 
above. 

The expression Eq. 1)5.12(1 for the spectrum can be 
summed analytically but yields an unsightly result which 
is not reproduced here. To compare with the experimen- 
tal results, we evaluate numerically the half- width at half 
maximum from S{uj). The results are plotted as a func- 
tion of probe power in Fig. 13 (full red line). The agree- 
ment with the experimental results (symbols) is good, 
especially given that there are no adjustable parameters 
apart from 72 which only sets the value of the dephasing 
at n = 0. In this figure we have included more experi- 
mental points, for higher powers, than presented in Fig. 
5 of the Letter. The experimental points presented here 
are obtained by fitting a Lorentzian (blue squares) and a 
Gaussian (red triangles) to the experimentally measured 
spectroscopic line (see Fig.EJ. We then keep the fit which 
has the smallest variance or both points if the variances 
are approximately the same (purple squares and trian- 
gles) . The error bars are the standard errors on the half- 
width half-max obtained from the fit. This approach to 
extracting the error bar is different from what was pre- 
sented in Fig. 5 of the Letter. In that case, the error 
bars represented the systematic difference between the 
Lorentzian and Gaussian fits and the points were the av- 
erage value. From Fig. O (this paper) we see that the first 
7 points fit best to a Lorentzian while the later points fit 
best to a Gaussian, except for the higher powers where 
the error in the fit is approximately the same. The pre- 
dicted crossover from Lorentzian to Gaussian is clearly 
seen. 

There are several potential sources of discrepancy be- 
tween the experimental and the theoretical results. One 



of them is the breakdown of the lowest order dispersive 
approximation. Using the same simple non-linear model 
as in section IIVI we plot in Fig. O the half- width at half 
maximum as a function of input power (green dashed 
line). The effect of this correction is to reduce the width 
of the spectroscopy peaks. 

The breakdown of the dispersive approximation can 
be seen by the dispersive result (red full line) over- 
estimating the width at high powers. However, while 
the simple non-linear model used here does correctly 
show saturation of the width at high powers, it is not a 
full treatment of the Jaynes-Cummings Hamiltonian and 
should not be considered too seriously. A complete inves- 
tigation of the behavior of the system at very large pho- 
ton numbers will require numerical investigation which 
is beyond the scope of this paper. 

A further possible source of discrepancy comes from 
the fact that a constant spectroscopy power was used. 
Since the effective coupling strength x is not constant 
with probe power, the effect of the spectroscopy power on 
the qubit will change with measurement probe power. As 
a result, spectroscopic power broadening will also de- 
pend on measurement probe power. This effect has been 
taken into account in the green dashed curve of Fig. |2| 
but only to the accuracy of our simple model. Finally, 
environmental noise due to two-level systems activated 
at large photon number could be an additional cause of 
discrepancy. 



B. Beyond the gaussian approximation 

To go beyond the gaussian approximation made in the 
last section, we solve the master equation Eq. (|2.6|l using 
the positive P- function method [31| ■ Following Ref. [33 , 
we first write the qubit-cavity density matrix as 

P = Pee\e){e\ + f)gg\g){g\+ peg\e){g\+ pge\g){e\, (5.13) 

where pij acts only in the cavity Hilbert space. As shown 
in appendix ^ this leads to four coupled differential 
equations for the operators pij . In the absence of qubit 
mixing (Ti processes), solving these coupled equations 
yields the time-evolved full density matrix 

pit) ^ c,e(0)|e)(e|® |«+(t))(a+(t)| 



Cagmg){g\^\a.it)){a.{t)\ 
Cegit)\e){g\(g>\a+it))(a^it)\ 
Cge{t)\g){e\®\a.{t))(a+it)l (5.14) 



where 



At) 



'eg 



(a_(i)|a+(0) 



(5.15) 



describes the decay of the qubit phase coherence. In the 
above expression, we have 



,(0)e- 



-i{u. 



-i'y2)tp-t2xJo a+{t')a'_it')dt' 



(5.16) 
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and age{t) = a*^g{t) with 



a+{t) = a\ + cxp[-(k/2 + ix + iAr)i][a+(0) - a; 

(5.17) 

-ie-ri / [n /2 + ix + i^r) ■ Moreover, we have 



where a'lj 
Q;_(t) = 



+ exp[-(K/2 -ix + iA^)i][a_(0) 



(5.18) 

with at_ — —ie-ci/ {k/2 — «x + lA^). In these expressions, 
q;j_ represent the steady-state value of the field (a) given 
that the qubit is either in its ground (— ) or excited (+) 
state. Recall that A,, is the detuning of the measurement 
beam from the bare cavity frequency. 

In the limit ^ 1 discussed previously, the decay of 
aeg{t) is given by 



it) 



where F,, 



Tn, = -2xlm[a;a! 



(0)exp[-(72 + r^)t] 



(n+ + n^)KX^ 
k2/4 + + A2 



(5.19) 



(5.20) 



is the generalized measurement-induced dephasing rate. 
In this expression, 
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FIG. 5: (Color online) Measurement-induced dephasing rate 
Fm as a function of the detuning between the bare res- 
onator frequency ujr and the measurement drive frequency 
u>d- The dephasing rate is divided by the measurement power 
elf since it only changes the overall scale and not the structure 
of Fm. The blue (dashed) line corresponds to the experimen- 
tal parameters given in section UTTl The red (full) line has the 
same parameters but a larger cavity pull x/Stt = 5 MHz. 



n± 



«;2/4+(A.±x)^ 



(5.21) 



is the stationary average number of photons in the cav- 
ity when the qubit is in the excited (+) or ground (— ) 
state. At small pulls, as is the case for the experimen- 
tal parameters quoted in section IIIII the measurement- 
induced dephasing rate is largest at A^ = (see Fig. El 
blue dashed line). At large pulls, cavity transmission 
decreases at A^ — and, as illustrated in Fig. [S] (red 
solid line), the maximum dephasing rate then occurs at 
A^ = ±-\/x^ ~ k2/4. As we increase x^ the information 
about the state of the qubit is conveyed more by the 
amplitude than the phase of the transmitted beam. 

From the above, we see explicitly that by introduc- 
ing a probe (Crf) we cause the coherence terms aeg{i) to 
exponentially decay at a rate F^ thereby leaving the sys- 
tem in a mixed state with perfect correlation between the 
eigenstates of Uz and the pointer states |a+) and \a-)- 
As a result, if the pointer states are well separated in 
phase space, we can regard the cavity as a meter which 
performs a von Neumann projective measurement of the 
qubit observable Uz- 

A measure of the distinguishability of the cavity states 
is £> = |q;+ — q;_|2 |31|. If this is large such that 
|(a_|a+)| = exp[— £)] is small, then the two pointer 
states are well separated and easily distinguishable. In 
the steady-state, D is given by 



2(nH 



kV4 + x' + A2' 
which is related to Fm in the following way 



(5.22) 



(5.23) 



That is, as the measurement becomes more projective in 
the dz basis, the qubit dephases faster. This is a clear 
example of measurement induced dephasing and of the 
fundamental limit which exists between acquiring infor- 
mation about a quantum system and dephasing of that 
system [s^. Note that, as shown appendixiBl can also 
be related to the measurement time. The present system 
is a factor of 4 away from the quantum limit. One factor 
of 2 comes from the fact that we are using a symmetric 
resonator and only looking at the transmission. Half of 
the information is lost in the unmeasured reflected sig- 
nal The other factor of two comes from the use of 
heterodyne rather than homodyne detection. As a result 
and as shown in this appendix, the quantum limit can be 
reached by using asymmetric resonators and homodyne 
detection of the transmitted field. 



1. The qubit absorption spectrum 

To evaluate the qubit 's spectrum, we first need to cal- 
culate the correlation function ((T_(t)(j+(0))s. Note that 
(as discussed in appendix ^) in calculating this partic- 
ular correlation function we do not need to assume that 
Ti is infinite. This is because it only depends on the 
off-diagonal coherences and these are not mixed by a Ti 
process. Using Eq. 15.2|l and the above results, the cor- 
relator can be shown to be {t > 0) 



(5.24) 
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with the initial condition Q!+(0) = a-(0) = aL. Using Eq. H5.25|l to order 
Eq. H5.16|l with the above initial condition yields 



X exphAe-(''/2+^'^+*^-'*]cxp[A], 



(5.25) 



where 



A 



B = 



(a! 



a , )a_ 



k/2 — «x — iAr 
k/2 + «x + iAr ' 
2xRe[a+al*] = x(^+ 



D. 



(5.26) 

n_)-x£'s. (5.27) 



From the above expression, we see that the time depen- 
dence of the correlation function is given by three terms. 
The first one involves Fm and is simply the Forentzian 
part of the measurement-induced dephasing spectrum we 
have seen before. The second is a frequency shift B which 
contains a negative term — x^s which gives rise to neg- 
ative frequency contributions in the spectrum. The last 
relevant term goes as Aexp[{K/2-|-i(x + Ar)}i] and gives 
rise to non-Lorentizian spectra. 

From the expression for the correlation function, it is 
simple to obtain the spectrum: 



1 °° 1 



TT ^ — ' 7 

3=0 



{-Aye^ 



Tj/2 - i{uj - ujj) 



(5.28) 



where Tj = 2(72 -I- F^) +jn and ujj 



-B + j{x + Ar 



The spectrum, as in the Gaussian approximation, can 
be written as a sum over different photon numbers j. 
In the limit of (A^. -I- x) much different from k/2, the 
spectrum is a sum of Forentzians with decay rate Fj/2. 
However unlike Eq. (|5.12() where each Forentzian is cen- 
tered at the ac-Stark shifted frequency, here each peak 
has its own frequency shift loj. As a result, the full the- 
ory predicts that the spectrum need not be symmetric 
whereas in the Gaussian theory only symmetric spectra 
are possible. Furthermore, in the limit that x is much 
larger than k, A —^ and the spectral weights become 
Poisson distributed with mean Ds- The peaks are sepa- 
rated by (x + Ar) and the first one is at the frequency 
Oa + B. That is, the average frequency, which is the 
ac-Stark shift, occurs at 



= B+{x + Ar)Ds = 2xn_ 



(5.29) 



Taking this limit further with x much larger than the 
widths Tj, the individual peaks will become distinguish- 
able. This is discussed further below. 

It is interesting to point out that in the limit of large 
n_ (or n_|_) the results obtained here and those obtained 
in the Gaussian approximation agree. This can be seen 
by expanding the exponent Aexp[{K/2 -I- i{x + Ar)}i] in 



(a_(t)a+(0)). 



(72 + icoa + i2xn-)t 



exp 

n ' . (5.30) 

i^(KV4+(A,. + x)^)|tp^ 



For Ar = 0, n_ = n+ = n and we recover Eq. 115.11(1 in 
the lab frame. As a result, in the large n limit, both the 
gaussian approximation and the above theory converge 
to give the same gaussian spectrum. It is only when n is 
small that the theories have different predictions. This 
will be discussed further in section lVBSI 

The half-width at half maximum of the spectrum ob- 
tained from the full expression Eq. (|^^ is plotted at 
Ar = as a function of n in Fig. |21 using the experi- 
mental parameters given in section Hill Since the experi- 
ment was done in the limit x < '*/2, the results obtained 
from Eq. (|5.28|) cannot be distinguished from those ob- 
tained from the Gaussian approximation Eq. H5.12|) . This 
is because in the small cavity pull limit Fm —f F,„ and 
therefore A — *■ 2F,n,/K, B 2nx- That is, to see the 
break down of the gaussian approximation at small n re- 
quires a larger cavity pull. 



2. Phase-gate 

As discussed above and in appendix ^ measurement 
causes dephasing of the qubit. This is a clear illustration 
of the Heisenberg type relation between rate of informa- 
tion gain and dephasing F^ [ssj . However, irradiation at 
the rf frequency does not have to induce dephasing of the 
qubit. Indeed, the qubit pulls the resonator frequency up 
or down causing a state dependent phase shift for pho- 
tons near the cavity frequency. But, in the low x limit, 
photons off resonant from the resonator have phase shifts 
nearly independent of the qubit state. These photons do 
not become entangled with the qubit, and hence do not 
cause dephasing. 

This can be understood more quantitatively in Fig. |3 
where it can be seen that the dephasing rate is significant 
only on a frequency range k around the pulled resonator 
frequency. In this Figure, we have fixed the input power 
and scanned A^. In Fig. El we rather keep the number 
of photons in the cavity fixed (n_ = 2) and scan A^. We 
see that at large detunings, the dephasing rate scales as 
A~^. As a result, off- resonant irradiation can produce 
large ac-Stark shifts (wac) with minimal dephasing of the 
qubit. This can be used as a single-bit phase gate for 
quantum computation. 

The observed asymmetry in the large x case (red solid 
line of Fig. Ejl is a result of the fact Fm depends on fi- 
and n+. By writing n+ as. 



^2/4 -t- (A. -x)^ 



(5.31) 



k2/4+(A. +x)2' 
we see that at fixed can be large for negative A^.. 
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3. Number splitting 
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FIG. 6: (Color online) Dephasing rate Fm as a function of 
the detuning Ar between the bare resonator frequency ujr 
and the measurement drive frequency ujd for fixed n_ = 2. a) 
Experimental parameters given in section UTTI b) Same as a) 
but with a larger cavity pull x/2n = 5 MHz. 



Thus the overall measurement induced dephasing will be 
large in that region. 

The quality factor for this single qubit gate, Q, can 
be defined as the coherent phase rotation that can be 
realized in the total dephasing time (Tj^^ + r„i)~^ 0. 
That is 



Q 



2(^2 



2kx 



(5.32) 



in the large pull limit and in the ideal situation where de- 
phasing is limited by photon shot noise. Moreover, sim- 
ilarly to Rabi oscillations that have been demonstrated 
experimentally this phase gate could be realized on 
a time scale which is much faster than 1/k since for off- 
resonant irradiation, the cavity is only virtually popu- 
lated. 

To show that the dephasing is minimal during the 
phase-gate, we have calculated using Eq. (|5.28() the 
linewidth of the qubit spectrum as a function of input 
power for the experimental parameters and a large pos- 
itive detuning Ar/27r 32 MHz (well away from the 
peak shown in Fig.|n|). This is shown in Fig. 13 as a blue 
dashed line. Here we see that in the dispersive model, 
there is no additional dephasing due to the off-resonant 
irradiation. That is, the predicted linewidth stays con- 
stant at 72 for the input powers plotted. At the critical 
photon number, ricnt = A^/4g^ the quality factor for 
the above experimental parameters is 17.3. This qual- 
ity factor can however be easily increased by optimizing 
the system parameters l^^l- For example at g/2TT = 100 
MHz and A/27r = 1000 MHz a quality factor of 157 can 
be reached. (This value is entirely limited by the current 
value of T2 ~ 500ns and not by the direct infidelity of 
the phase gate.) An advantage of this rf approach over 
a dc pulse of the flux or gate charge is that the logical 
operation can be realized while biased at the sweet spot. 
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FIG. 7: (Color online) Spectrum S{u) as a function of detun- 
ing At/k, at fixed cavity pull x/ hi = 100. The dephasing rate 
was set to the conservative value of 72 = 7.6k. The average 
photon number in the cavity was fixed to 7i_ = 2 in panel 
a) (Ar > 0) and we choose n+ = 2 in panel b) (A,- < 0). 
This implies that the measurement beam power changes with 
detuning. Insert: Spectrum at Ar = X = 100k. At this de- 
tuning, the peaks are split by 2x- More generally, they are 
split by X + A. Large detuning yields large splitting but, as 
shown in Fig. |K1 this can be at the expense of small measur- 
able phase shift in the transmitted field. 



As explained in the previous sections, the gaussian 
model and the P-function approach agree in the small 
pull limit X ^ '^/S- In the large pull case, the predicted 
behavior is however substantially different. Indeed, when 
the Lorcntzians in Eq. H5.28|l are separated in frequency 
by more than their width, the spectrum S'(w) will be 
split into many peaks with each peak corresponding to 
a different photon number in the cavity. Number split- 
ting was also predicted by Dykman and Krivoglaz for a 
different situation, namely an undriven cavity coupled 
to a thermal bath [s^- The number splitting for our 
case (driven cavity at zero temperature) is illustrated in 
Figures [T] and |S1 In Fig. [7| the spectrum is shown as a 
function of frequency and of the detuning for a fixed 
x/k, = 100. In Fig. [7^) (Ar > 0), n_ is fixed and equal 
to 2 whereas in Fig.[7|D) (A^ < 0), n+ = 2. The choice of 
fixing either or n_ was made so that there is always a 
small number of photons in the cavity independent of the 
state of the qubit. For example, for n_ fixed, the number 
of photons in the cavity when the qubit is excited, n_|_, 
would be small for A^ > 0, equal to n_ at A^ = and 
very large at A^ = — x- The inset shows a cross-section 
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Frequency, (co-coa-x)/27i [MHz] 

FIG. 8: (Color online) Spectrum S{u!) as a function of the 
cavity pull x/^- The detuning is Ar — x such that the cavity 
is always driven at the dressed frequency uj^ — X- Other pa- 
rameter values are the same as in Fig. 13 Inset: Spectrum at 
xl = 20. At this experimentally realistic value of the cavity 
pull, the number splitting should be resolvable. 



at Ar/K — 100 (Ar — X). For this value of A^/k, the 
peaks are very well separated and the integrated area un- 
der each peak obeys Poisson statistics. It is interesting 
to stress that only at A^ = x does the spectrum have a 
simple Poisson distribution corresponding to a coherent 
state with average photon number = 2. At A^. < x 
there arc more peaks than expected for a coherent state 
of this amplitude. Furthermore, at A^ = — x the spec- 
trum is single peaked. 

In Fig. |H| the spectrum is shown as a function of the 
cavity pull x/k for a detuning of A^ ~ X; such that the 
cavity is always driven at the pulled frequency lOr — X 
corresponding to the qubit in the ground state. In these 
plots, we have taken 72 = 7.6k. Assuming an exper- 
imentally realistic value of /t/27r ~ 100 KHz |23|, this 
corresponds to a conservative T2 = 200 ns 0. As seen 
on Fig |S1 for these parameter values, the peaks should 
be resolvable experimentally starting around xl ~ 20 
(insert in Fig. IS)). Achieving x/k ~ 20 in the dispersive 
limit (5/A < 0.1) requires g/27r ~ 20 MHz. This value 
of g was already realized experimentally and therefore 
the experimental observation of number splitting seems 
feasible. 

The behaviour described above can be understood sim- 
ply as ringing of a high-Q resonator when its resonance 
frequency is suddenly changed when the qubit changes 
state. Equivalcntly we can think of this as a Raman pro- 
cess in which drive photons in the cavity at the time of 
the transition are lifted up to the final cavity frequency. 



As discussed previously, the calculation of the spectrum 
assumes that the qubit is initially in the ground state 
with the measurement beam turned on at a frequency lo^ 
detuned by A^ from the bare resonator frequency w,.. In 
the calculation of the correlation function ((7_ (i)(T+(0))s, 
the qubit is flipped to the excited state at time i = and 
the overlap with the ground state is calculated at time 
t. When the qubit is flipped, the dressed resonator fre- 
quency is suddenly changed from — X to tj^ + X- De- 
pending on the frequency of the measurement drive and 
the quality factor of the cavity, this sudden change will 
cause ringing in the cavity. This is illustrated in Fig. |3 
where the distance D{€) = \aj^{t) — a_(t)p is plotted as 
a function of time for two values of xl For the mod- 
erate value of x/'* = 5 (red, full line), the distance is 
seen to undergo large oscillations before settling to the 
steady-state value D^. For the low x/k ratio of 0.1 (blue, 
dashed line) there is no ringing due to the abrupt change 
of dressed cavity frequency and the distance simply rises 
to A- 

The ringing is also shown in the inset of Fig. where 
the real and imaginary part of the cavity field oij^(t) are 
plotted as a function of time, again for xl ^ = 0.1 and 
5. In the low Q case, as the qubit flips, the cavity field 
settles to its new steady state value without large excur- 
sions in the held amplitude, and therefore large changes 
in photon number. Only a few different photon num- 
bers contribute and the corresponding spectrum is single 
peaked as expected. In the high Q case, the field ampli- 
tudes performs many cycles before settling to the steady- 
state value. The cavity therefore probes a large range of 
photon numbers and the spectrum shows multiple peaks. 
In the time domain we can see that the qubit correlator, 
Eq. H5.25(l . has period recurrences provided x 2> k. It is 
these recurrences in time which give peaks in the spec- 
trum. These time domain recurrences were observed in 
the resonant regime (A = A,. = 0) with Rydberg atoms 
in Ref. |2|. 

In the case where A^ = — x {'^d ~ +x)i the cavity is 
driven at the dressed cavity frequency corresponding to 
the qubit in the excited state. In this situation, flipping 
the qubit does not produce any inelastic Raman scatter- 
ing (ringing) since the photons are already at the final 
cavity frequency. This is seen in Fig. at A^/k — —100 
where the spectrum is single peaked. 

Remarkably when the detuning is such that A^ < X; 
and as can be seen in Fig. [7| the spectrum has peaks 
at negative frequencies (i.e. below uja) as well as more 
peaks than expected. In particular, in the limit that 
X 3> K and at A,. = 0, the frequencies start at uja — 
2f7,_x with peak separation x and the spectral weights 
have a poisson distribution with mean 4fi and not n. 
To understand the presence of these peaks, we move the 
dispersive Hamiltonian Eq. H2.5|l to the frame defined by 
the unitary operator 

U ^tl+D[a+]+Il^D[a^]. (5.33) 
Here D[a\ is the displacement operator for the cavity 
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FIG. 9: (Color online) Distance D{t) = \a+{t) - a^{t)f as 
a function of time for Ar — x Erf ~ 2k. The red (full) 
line corresponds to x/^ = 5 and the blue (dashed) line to 
x/k = 0.1. Other parameter values are the same as in Fig.|3 
Inset: Real and imaginary part of the cavity field Q+(t) as a 
function of time. The arrows indicate direction of time. 



defined by 

D[a] = exp[Q;a^ — a*a], (5.34) 

and n± are the projectors for the excited and ground 
state of the qubit. That is, we move to a frame that 
takes both the pointer states out of the picture. In this 
frame, Eq. 12.5|l becomes 

H ^ -[(^a - '2xn]^z + fixa^aa;:. (5.35) 

Here we have considered the situation where ~ and 
neglected the last term of Eq. 1)2.5(1 which only leads to a 
small shift of the qubit transition frequency in the present 
situation. This Hamiltonian is pictorially represented in 
Fig. ^1 We immediately see that, in this frame, the 
qubit transition frequency is reduced by 2xn from the 
lamb shifted frequency Wa- Moreover, when the qubit 
is in the ground state, the Hamiltonian corresponds to a 
shifted and inverted harmonic oscillator (LHS of Fig. ^J. 
On the other hand, when the qubit is in the excited state, 
the harmonic oscillator is shifted but not inverted (RHS 
of Fig. lTUI) . Starting with the qubit initially in the ground 
state and the field in a coherent state of amplitude a_ 
corresponds, in the frame defined by Eq. H5.33|l . to a qubit 
in the ground state and the vacuum state of the oscillator. 
In this frame, flipping the qubit at time t = corresponds 
to applying the operator (for = 0) 

U'la+U oc&+D[a^-a+]=a+D[2Vn\. (5.36) 

The result is to both flip the qubit and to displace the 
oscillator to a coherent state of mean photon number 
An. This is exactly the observed structure in the qubit 
power spectrum and each of the observed peaks corre- 
sponds to one of the possible transitions between these 
two oscillators. We note that the above is similar to the 
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FIG. 10: (Color online) Pictorial representation of the disper- 
sive Hamiltonian in a frame defined by Eq. 15.331 at Ar — 0. 
When calculating the correlator {a-{t)a+{0)) , the qubit is 
flipped at time t = 0. In the displaced frame, this corre- 
sponds to both flipping the qubit and displacing the oscillator 
from the vacuum state to a coherent state with mean photon 
number 4n. This distribution is shown on top of the harmonic 
oscillator corresponding to the qubit in its excited state. This 
simple picture explains the observed peaks in the qubit power 
spectrum and the presence of negative shifts. 



MoUow triplet, where transitions both above and below 
the atomic transition frequency are possible due to dress- 
ing of the atomic levels by the presence of a strong pump 
drive [soj . 

Observation of number splitting would constitute a 
simple test of number quantization of the field inside the 
resonator in the dispersive regime. In the resonant regime 
(A = Ar ~ 0), number quantization has been verified in 
cavity QED using Rydberg atoms This was done 

by looking at the Fourier components in the probabil- 
ity to find the atom in the excited state as a function of 
time. In the dispersive regime (A > g, k) the recurrence 
time for the cavity field was too long to be observable 
[2^ . In the circuit QED system it should be possible to 
reach the strong dispersive limit x > '^i l/?2 where it is 
possible to probe number quantization. In the dispersive 
regime the qubit spectrum acts as a probe of the cav- 
ity field. When the rate at which information about the 
cavity state is passed to the qubit faster than the rate at 
which the cavity state changes significantly due to damp- 
ing and qubit dephasing, it is possible to learn about 
the statistics of the field from the qubit spectrum. Note 
the above predictions are only valid in the limit where 
X^/A < K. When this in not the case the higher order 
effects in the Jaynes-Cummings Hamiltonian [Eq. (|2.3|l ] 
will become important and will lead to non-Poissonnian 
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statistics of the resonator field. This will be discussed fur- 
ther elsewhere. Finally, we note an interesting proposal 
by Brune et al. [2^ l3nl | and a recent experiment |^ to 
prepare a Fock state of the cavity field (number squeez- 
ing) containing a single photon by monitoring the state 
of a continuous beam of atoms sent through the cavity. 



VI. CONCLUSION 

We have found that due to the ac-Stark shift, quan- 
tum noise in the photon number populating the resonator 
leads to well characterized measurement-induced dephas- 
ing. A simple model based on a gaussian approximation 
for the phase noise was presented, as well as a more gen- 
eral model based on the positive P-function. For the 
experimental parameters given in the Letter, both mod- 
els yield the same quantitative results which are in very 
good agreement with the experiment. We emphasize that 
the only adjustable parameter in the theory is the intrin- 
sic qubit dephasing rate whose only effect is to give a 
constant offset to the predicted linewidth. 

In the strong dispersive regime, where the cavity pull x 
is much bigger than the cavity field decay rate k/2, the P- 
function approach predicts a splitting of the qubit spec- 
trum due to the discrete quantum nature of the field pop- 
ulating the cavity. Observation of this prediction would 
be a confirmation of the quantized nature of the resonator 
field. This strong dispersive coupling regime should be 
readily achievable with realistic circuit QED parameters. 
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APPENDIX A: POSITIVE-P REPRESENTATION 

In this appendix, we show how to solve the master 
equation Eq. H2.6|l in the presence of a measurement 
drive but take the spectroscopy drive Eg = 0. Substi- 
tuting Eq. 15.13|l into Eq. (|2.6|l yields the following four 



coupled differential equations 

Pee = n'D[a]pee - JlPee - t^rf [S + a\ Pee] - «X[a^a, Pee] 
- iAr[a^ a, Pee], 

(Al) 

Pgg = K'D[d]Pgg + 7l/5ee " i^rf + a\pgg] + IXia^CL, /5gg] 

— iAr[a}a, Pgg], 

(A2) 

Peg =K'D[a]peg -l2Peg - ieyl[a + o) , Peg] ~ ix{a^ a, Peg} 

- iAr[a^a, Peg] - iuJaPeg, 

(A3) 

Pge = KT>[a]pge - 72Pge - «£rf [S + a\ Pge] + ix{a^ a, pge} 

- iAr[a^a, Pge] + iCbaPge- 

(A4) 

Solving these four differential equations would yield a 
complete solution. However because of the coupling in- 
troduced by 71, this is not possible analytically for all 
possible observables. However for the particular case of 
computing the dephasing rate, we can (without error) set 
71 = in the equations for pee and pgg while keeping the 
contribution of relaxation to dephasing in the equations 
for the off-diagonal components. 

To solve Eqs. (|Il|l - (|X4jl we express the density matrix 
under the positive P-representation |3ll |: 

p,, = d^/3^^P.,(a,/?). (A5) 

Using this expression in Eqs. (|Al|l - (jA4|l and the identi- 
ties 



a\a) = a|a), 

a)\a) = (a„+a72)|a), 

(/3*|at = /3(/3*|, 

(/3*|a = {dp + P/2){P*l 



(A6) 
(A7) 
(A8) 
(A9) 



gives four coupled differential equations for the 'proba- 
bility densities' Pij: 



Pee = da [(i^rf + iX^ + lA^a + Ka/2)Pee] 

+ dp[{-ie,f - ixP - iArl3 + Kp/2)Pee] 

Pgg = da [(j^rf ~ ^X" + lArU + Ka/2)Pgg] 

+ df3[{-ie,f + 1x13 - iAr/3 + KP/2)Pgg] 



(AlO) 
(All) 



Peg = da[{ierf + ix^ + lArtt + Ka/2)Peg] 

+ dp [(-i£rf + ixP - I Ar 13 + KP/2)Peg] (A12) 
- i2xa(3Peg - l2Peg ~ il^aPeg 

Pge = da[{ie.rf - IXO. + lAr-a + Ka/2)Pge] 

+ dpii^isrt - ixl3 - iArf3 + Kl3/2)Pge] (A13) 

-I- i2xaPPge - 72-Pge + ioJaPge- 
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To obtain these expressions, we have assumed that 
Pij{oo,oo) = as is usual [25l |. 

These equations can be solved simply by making the 
ansatzen 



Pee ^S^'^C 



Peg -fleg W^''^ [« - a+{t)]5^^^ [/? - a*L{t)], 
Pge [« - a- [f3 - aUt)] 



(A14) 



and substitute these into each equation. This results in 



leg 



a, 



ge 



= -i£rf-«(A,. + x-«K/2)a+ (A15) 

= — i£rf — i — X ~ *^/2) «- (A16) 

^ - i {iva - ij2) aeg - i2xa+a*_aeg (A17) 

=i (a)a + 172) flge + i2xa_a+age- (A18) 



Solving these simple differential equations completely 
solves the master equation Eq. (|2.6() in the absence 
of mixing due to Ti effects. The solution is given in 
Eqs. - fTTFIl . 



From this quantum trajectory the rate at which in- 
formation is obtained about (a) is Krj. To convert this 
to a rate of information gain about (a^) we define the 
measurement observable for a time r as 



/(r) = 



1 



Re[J{t)e-'^]dt, 



(B6) 



where determines the quadrature in which the informa- 
tion about the qubit is encoded. This can be determined 
by 



tan 6 — 



Im[a^ 



Refa^ 



(B7) 



For example, if Ay, = and = 0, the information about 
the qubit is only encoded into the real part of (a) . From 
this observable, the mean and variance is 



I{t) ~ K77Re[(a)e 
A|/(r)| = 



[/(r) - /(r) 



(B8) 
(B9) 



APPENDIX B: MEASUREMENT TIME 

In this appendix we show how one can calculate the 
measurement time for this system and show how it re- 
lates to the quantum limit j35| . To do this we need to 
describe how we are measuring the pointer states (i.e. 
how the information is processed). In this experiment, 
this is done by using heterodyne detection of the signal 
that is transmitted from the cavity. In other words, the 
full quantum trajectory for the system is |40ll4l]| 

dpj{t) = dtCpj{t) + dm[{r{t) - Kri{d)))a\pj{t), (Bl) 

where Ti.[{J* {t) — K,rj{a^))d]pj{t) is the superoperator rep- 
resenting the non-linear effects of the continuous moni- 
toring and is defined by 



n[A]p = Ap + pit - (i + it)p. 



(B2) 



The measurement record (heterodyne signal) is given by 



J{t) = KTj(d) + ^C{t), 



(B3) 



where (^(t) is a complex gaussian white noise term, which 
is formally defined as 



mmt')] = mt)] = o 
mt)cit')] = s{t-t') 



(B4) 
(B5) 



where E denotes an ensemble average and 77 = 1/[2(A^ -|- 
1)] is the inefficiency of the measurement. Here N is the 
dark noise and the extra factor of 1/2 is due to the fact 
that information leaks out of the cavity in both directions 
and we only monitor transmission 



vity ] 

m 



That is, if we were to measure the system for a time t 
many times we are confident that to one standard devi- 
ation the value of I is K77Rc[(a)e~*"^] ± y/ Kr]/2T. 

If r is much shorter than I/7 then we can approximate 
I{t) with I±{t) = KriKc[a-^e^^'^], where the ± subscript 
refers to the state of the qubit. To be able to distinguish 
between (ct^) = ±1, we require A{az) < 1 and thus 



Al/I < 



U 



The equality defines the measurement time 
the above, this can be rewritten as 



1 



(BIO) 
Using 

(Bll) 



That is, even for perfect detection efficiency 77 = 1, this 
approach is a factor of two away from the quantum limit 
^measrm = 1/2 ^| . This is becausc even though we are 
selecting the correct quadrature in which the informa- 
tion about the qubit is stored, [using the classical pro- 
cessmg defined in Eq. IpTjl ] we are still measuring the 
other quadrature as we are performing heterodyne detec- 
tion. It is well known that heterodyne detection measures 
both the (/) and 4> -\- 7r/2 quadrature with 1/2 efficiency 
[4^ |4^ . Thus if we change the detection scheme to ho- 
modyne detection of the quadrature we can reach the 
quantum limit. That is, to reach the quantum limit we 
require 77 = 1 which means we need asymmetric cavities 
and no dark noise as well as a detection scheme which 
extracts only information about cTz- Note if we did not 
perform any classical processing on the heterodyne sig- 
nal J{t) then we would be a factor of four away from the 
quantum limit. 
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